Ground-plane screening of Coulomb interactions in two-dimensional systems: How 
effectively can one two-dimensional system screen interactions in another? 
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The use of a nearby metallic ground-plane to limit the range of the Coulomb interactions between 
carriers is a useful approach in studying the physics of two-dimensional (2D) systems. This approach 
has been used to study Wigner crystallization of electrons on the surface of liquid helium, and most 
recently, the insulating and metallic states of semiconductor-based two-dimensional systems. In this 
paper, we perform calculations of the screening effect of one 2D system on another and show that 
a 2D system is at least as effective as a metal in screening Coulomb interactions. We also show 
that the recent observation of the reduced effect of the ground-plane when the 2D system is in the 
metallic regime is due to intralayer screening. 

PACS numbers: 71.30.+h, 71.10.-w, 71.45.Gm 
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I. INTRODUCTION 

In a two-dimensional electron system (2DES), strong 
Coulomb interactions between electrons can lead to ex- 
otic phenomena such as the Wigner crystal state,— i^i^ the 
fractional quantum Hall effect ^''^ and the anomalous 2D 
metallic statej^iii^ One route to studying the role played 
by Coulomb interactions is to limit their length-scale us- 
ing a metallic ground-plane located close to the 2DESi^iiS 
This approach was first used in studies of the melting of 
the Wigner crystal state formed in electrons on a liquid 
He surfaceiiiii^ More recently, it has been used to study 
the role of Coulomb interactions in the insulatingi^ and 
metallioi^ regimes of a 2D hole system (2DHS) formed 
in an AlGaAs/GaAs heterostructure. 

Whereas the study of Coulomb interactions in the in- 
sulating regimei^ was achieved quite straightforwardly 
using a metal surface gate separated from the 2DHS by 
~ 500 nm (sec Fig. [TJa)), the corresponding study in 
the metallic regime could not be achieved in this way. 
This is because the higher hole density p in the metallic 
regime requires that the distance d between the 2DHS 
and ground-plane be comparable to the carrier spacing 
{d ^ 2(7rp)^^/^ ^ 50 nm) to achieve effective screening, 
and at the same time that the 2DHS be deep enough 
in the heterostructure (> 100 nm) to achieve a mobility 
sufficient to observe the metallic behavior. To overcome 
this challenge, a double quantum well heterostructure 
was used (Sec Fig. [IJb)) such that the 2DHS formed 
in the upper quantum well (screening layer) served as 
the ground-plane for the lower quantum well (transport 
layer), enabling the measurement of a ~ 340 nm deep, 
high mobility 2DHS separated by only 50 nm from a 
ground-plane.— 

In considering experiments on screening in double 
quantum well systems, a natural question to ask is 
whether a 2D system is as effective as a metal gate when 
used as a ground-plane to screen Coulomb interactions 
between carriers in a nearby 2D system. This is impor- 



tant given that the screening charge in a 2D system is 
restricted to two dimensions and the density of states is 
several orders of magnitude smaller than in a metal film. 
In this paper, we perform calculations of the screening ef- 
fect of a ground-plane on a 2D system for two cases: The 
first where the ground-plane is a metal and the second 
where it is a 2D system. We begin using the Thomas- 
Fermi approximation in the absence of intralayer screen- 
ing in the transport layer to show that a 2D system is at 
least as effective as a metal gate as a ground-plane for 
the experiment in Rcf. 14. We also compare the experi- 
ments in the insulating^ and metallioi^ regimes of a 2D 
hole system in the Thomas-Fermi approximation, to ex- 
plain why the ground-plane has less effect in the metallic 
regime compared to the insulating regime. Finally, since 
the experiment by Ho et al. was performed at Tg > 1, 
where the Thomas-Fermi approximation begins to break 
down, we extend our model to account for exchange and 
finite thickness effects to see how these affect the conclu- 
sions from the Thomas-Fermi model. 

The paper is structured as follows. In Section II we 
derive the dielectric functions for screening of a 2D hole 
system by a metal gate and another nearby 2D hole sys- 
tem. In Section HI, we compare the various dielectric 
functions numerically and discuss their implications for 
the ground-plane screening experiments of 2D systems 
in the insulatingi^ and metallioii. regimes. Conclusions 
will be presented in Section IV. For readers unfamil- 
iar with the intricacies of screening in 2D systems, we 
give a brief introduction to the screening theory for a 
single 2D system in Appendix A to aid them in under- 
standing the theory developed in Sect. II. In Appendix 
B, we compare our model accounting for exchange and 
finite thickness effects to related works on many-body 
physics in double quantum well structures. In Appendix 
C we show how our analysis of ground plane screening 
can be related to previous work on bilayer 2D systems 
that have studied the compressibility. In contrast to ex- 
periments on ground-plane screening of the interactions 
within a 2D system^Hii^ previous experiments on bilayer 
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FIG. 1: 



Schematics of the ground-plane screening experiments 
recently performed by (a) Huang et alM- and (b) Ho et ali^. 
In (b), there are two possible ground-plane configurations. 
In the first, the gate is grounded and the 2D system acts as 

the ground-plane. In the second, the gate is biased to 
deplete the upper 2D system, and the gate then acts as the 
ground-plane instead. This allows the distance between the 
transport layer and the ground-plane to be varied in situ - 
For more details, see Ref. 14. 
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FIG. 2: 

Schematics showing the two systems considered in this 
paper. The transport layer is screened by (a) a metal surface 
gate and (b) a second 2D system. In both cases the 
screening layer is separated by a distance d from the 
transport layer, and the transport (1) and screening (2) 
layers have independent potentials and charge densities p. 



II. THE SCREENING OF ONE 2D LAYER BY 
ANOTHER 



2D systems have examined the compressibility of a 2D 
system by studying penetration of an electric field ap- 
plied perpendicular to the 2D planed. We show that 
although these two screening configurations have quite 
distinct geometries, our model is consistent with widely 
used analysis^ describing the penetration of a perpen- 
dicular electric field through a 2D system. 

In the calculations that follow, we use linear screening 
theory and the static dielectric function approximation 
(i.e., uj 0). Unless otherwise specified, we assume for 
convenience that the 2D systems contain holes (electron 
results can be obtained with appropriate corrections for 
charge and mass) to facilitate direct connection with re- 
cent experimental results in AlGaAs/GaAs heterostruc- 
turesJ^ii^ We also assume that tunneling between the 
two quantum wells is negligible and ignore any Coulomb 
drag effects (i.e., interlayer exchange and correlations). 



We now begin considering the screening effect of a 
nearby ground-plane on a 2D system (transport layer) 
for two different configurations. In the first, the ground- 
plane (i.e., screening layer) is a metal surface gate (see 
Fig.[2ja)) and in the second, the ground-plane is another 
2D system (see Fig. [H^b)). In both cases the transport 
and screening layers are separated by a distance d. 

If we consider some positive external test charge p^^* 
added to the transport layer, this leads to induced charge 
in both the transport layer p™'* (as in Appendix A) and 
in the screening layer p™''- Note however that no external 
charge is added to the screening layer, so pj^* = and 
P2 ~ pT"^, whereas pi = pf^* + p™''. 

How we deal with the induced charge in the screening 
layer differs in the two cases. In both cases, we consider 
the transport layer and a 2nd layer of charge a distance D 
above it; each layer having a potential and charge density 
of 4>i{q), pi{q) and 4>2{q), P2{q); respectively. For a metal 
surface gate, we can use the standard image charge ap- 
proach which involves considering the induced charge 
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in the screening layer as a 2D layer of negative 'image' 
charges located a distance D = 2d away from the trans- 
port layer. This results in p™'' = — pi for a metal gate. 
The image charge approach assumes that the ground- 
plane is a perfect metal. This assumption is relatively 
well satisfied by a typical metal surface gate (Au gate 
~ 150 nm thick) but not by a 2D system. Thus when 
a 2D system is used as the screening layer, we cannot 
assume an added induced negative image charge as we 
can for the metal. Instead, we account for screening by 
a 2D system by directly calculating its induced charge, 
which is located in the screening layer (i.e., at a distance 
D = d). This results in p™'* = for a 2D screen- 

ing layer, where x'^ is the polarizability, which describes 
how much p™'^ is produced in response to the addition of 
pext^ Note that we are only considering the case where 
the both 2D layers are of the same type of charge (e.g., 
a bilayer 2DHS). 

Charge in one layer leads to a potential in the other 
via the interlayer Coulomb interactioni^^ 
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V{q) 



(1) 



where V{q) 



2^ is the intralayer Coulomb interac- 
tion, D is the distance between the two layers, and 
is the Fourier transform. If the screening layer is a metal 
D = 2d, and Z? = d if it is a 2D system. The resulting 
potential in the transport layer then becomes: 



M^)^V{q)p, + U{q)pT''in) 



(2) 



We discuss how to obtain P2"'(ct) section II. A. The ef- 
fectiveness of the screening is obtained from the modified 
dielectric function e{q,d), which we define as the inverse 
of the ratio of the screened potential to the unscreened 
potential for the transport layer: 



Thomas- Fermi (TF) approximation^^ At first we will ig- 
nore (Sect. IIA) and later include (Sect. IIB) intralayer 
screening in the transport layer in the calculations. Fi- 
nally, in Sect. IIC, we will extend the model for a 2D 
screening layer to account for its behavior at lower den- 
sities, to confirm that our conclusions from the simpler 
calculations arc robust. 



A. No Intralayer Screening in the Transport Layer 

We begin by considering the case where there is no in- 
tralayer screening in the transport layer. This is useful 
because it allows a straightforward comparison of the ef- 
fectiveness of the 2DHS as a ground-plane, without the 
obscuring effect of intralayer screening. To do this cal- 
culation, we set p™"* = 0, such that pi = pf^*. In other 
words, there is only external charge in the transport layer 
and only induced charge in the screening layer. 

Considering the metal gate first, we have p™'^{q) ~ 
~Pi{l) = ~Pi^'(q) from the method of images. If 
we combine the two results above for pi and p™'' with 
Eqn. [21 we obtain: 



M^) = iV{q) - Uiq))priq) 



(4) 



After using Eqn. [T] to eliminate U{q), Eqn. [3] then gives 
the dielectric function for the metal gate: 
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Cmctal.ns 



(5) 



where the additional subscript ns denotes that intralayer 
screening has been ignored. 

For a 2D screening layer, the dielectric function is ob- 
tained self-consistently through the RPA as follows. The 
induced charge is related to the screening layer potential 
by: 
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(3) 



The dielectric function for the transport layer can be 
obtained for three possible configurations: no screening 
layer (i.e., just a single 2D system), a 2D screening layer 
and a metal screening layer, which we denote as esingio('i'), 
e2Y){q,d) and emetai(9, f^), respectively. esingio('?) is the 
conventional dielectric function for a single 2D system, 
and can be recovered by taking the limit d ^ oo for 
either e2D('7,ci) or CmotaK'?, c^), as we show in Section II 
Part C. If one wanted to take transport theory for a single 
layer, and consider the effects of ground plane screening, 
it is only necessary to replace esingie(9) with t2T>{q,d) or 
emctai(9, d) to take into account the screening effect of a 
2D or metal ground plane, respectively. 

We will now obtain the various e using an approach in- 
volving the Random Phase Approximation (RPA)ii and 



02(q) = C/(g)pi + F(q)pr'^(q) 
P2"'(q) = X^(9)02(q) 



(6) 
(7) 



where x% is the polarizability of the screening layer, nor- 
mally given by the 2D Lindhard functioni^^ When this 
is combined with Eqn. [51 knowing that pi = pf^*, we 
obtain: 



ind _ xliqWiq) -qd 



(8) 



This result is substituted into Eqn. [2l and using Eqns. [T] 
and [3] gives: 



e2D,ns(g, d) 



= 1 



xl{q)V{q) 

^~xl{q)v{q)' 



-2qd 



(9) 
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To simplify this expression, we use the Thomas-Fermi 



,2 dn 



approximation X2il) — dJI 
dynamic density of states of the 2D system,— 



where ^ 



is the thermo- 
to give: 



e2D,ns('7, 



q + q 



TF 



-2qd 



TF 



(10) 



where the Thomas-Fermi wavevector q — t^t. 

Note that if wc take the 2D screening layer to the metal- 
lic limit, in other words, wc give it an infinite density 
of states, which corresponds to q^^ oo, then Eqn. [TOl 
reduces to Eqn.O as one would expect. 



With Intralayer Screening in the Transport 
Layer 



the transport layer. To approach this problem, we again 
place an external charge density pf^* in the transport 
layer, but now we have induced charge in both the trans- 
port p™'' and screening p™'' layers. Additionally, we la- 
bel the polarization Xi{q) ^^^d Thomas-Fermi wavenum- 
ber qf^ where i = 1 or 2 corresponding to the transport 
and screening layers, respectively. The derivation pro- 
ceeds as before, but with the addition of the induced 
charge density Pi"''^(q) = Xi(9)0i(q) i^i the transport 
layer. The results obtained are: 
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1 



-2qd 



(q.d) 



i-^(9)x?(g)(i-e-2?'i) 

1 — e^^'^^ 



1 + ^(1 

a ^ 



-2qd^ 



(11) 



We now consider the case where there is intralayer 
screening (i.e., finite polarizability and induced charge) in 



for the metal gate, where the added subscript s denotes 
that intralayer screening has been included, and: 



1 - VxVi^ - e-'^'''^] 



1 + 



^-2qd] 



e2D,.(g,d) [1 - Vx\][l - VxD] - V^xlxle-^'^" 



(1 + V)(i + V) 



(12) 



when the screening layer is a 2D system. 

We can check the consistency of these equations with 
those in Sect. IIA in three ways. Firstly, by setting 
qi^ = 0, which corresponds to no screening or induced 
charge in the transport layer, Egns. [TT] and [T^ reduce to 
Eqns. O and [TUl respectively. Secondly, if we set q^^ = 
instead, which corresponds to no screening or induced 
charge in the screening layer, then s Eqn. [T^] re- 
duces to e^,j which is given in the Thomas-Fermi ap- 
proximation byi^ 

1 TF 

= + (13) 
esingio(q) q 

Finally, if we set qj^ ^ oo to take the 2D screening 
layer to the metallic limit, then ^ in Eqn. I12l reduces 

to emctai,s in Eqn.ini 

C. More Accurate Calculations for 2D Systems at 
Lower Densities 

Following our relatively simple treatment of ground- 
plane screening above, it is now interesting to ask how the 
results of our calculations change if we extend our model 
to account for two phenomena ignored in our Thomas- 
Fermi model: exchange effects at low densities, and the 
finite thickness of the screening and transport layers. 



The Thomas-Fermi approximation works well when 
the interaction parameter = {0'*b^/t^)~^ ^ 1, where 
a*g ~ AiTef^ / m* is the effective Bohr radius. However, 
it is not as accurate for 2D systems at lower densities, 
such as those used in our experiments^, where the inter- 
action parameters for the screening and transport layers 
were ^ 10 and 10.2 < rg < 14.3, respectively. At such 
low densities, it is essential to include the effects of ex- 
change, and a better approximation involves the addition 
of the local field corrections^, which we will use in this 
section. 

One might question if linear screening theory is valid 
at the large values considered here, or whether it is 
necessary to include non-linear screening effects. Recent 
compressibility measurements of 2D systems across the 
apparent metal-insulator transition^ have shown that 
the divergence of the inverse compressibility at low den- 
sities is consistent with non- linear screening theories j^i^ 
which take into account a percolation transition in the 2D 
system. This consistency suggested that at sufficiently 
low densities, non-linear theories were required to ade- 
quately describe screening. However, although we are 
interested in fairly high Ts, our analysis relates to mea- 
surements performed on the metallic side of the apparent 
metal-insulator transition, where the effect of inhomo- 
gcneities is minimal. Indeed, Ref. 25 shows a comparison 
of the non-linear screening theorie a^'^'^^ with a standard 
uniform screening theory. These theories all yield the 
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same results in our density range of interest, confirming 
that it is safe to treat the 2D system as homogenous and 
use hnear screening in our calculations. 

Returning to the inclusion of the local field correction 
to our calculations; when considering the case of two 2D 
layers, we will use the single layer local field factor G(q) 
to account for intralayer exchange effects, and for sim- 
plicity, ignore any corresponding intcrlayer effects. In 
this work we will use the Hubbard approximation for 
G(q) (see Egn. IA6l in Appendix A). This leads to: 



confined to 20 nm wide quantum wells in Rcf. 14. To 
do this, we introduce a form factor F{q) that modi- 
fies the bare Coulomb interaction such that V{q) 
V{q)F{q)^ The form factor is defined as F{q) = 

/ J \ip{z)\'^\tp{z')\'^e~''^^~^ ^dzdz' , where iIj{z) is the wave- 
function of an electron/hole in the direction perpendic- 
ular to the plane of the quantum welli^ Assuming an 
infinite-square potential for the quantum well, we ob- 
taini^i^ 



x^iq) 



X°(9) 



l-Viq)xnm-G^{q)] 



(14) 



where Gi{q) is the local field factor for layer i = 
(1,2), respectively. The calculations proceed as before 
in Sect. IIB, except that where we consider intralayer 
screening in the transport layer, we have: 



1 ro Stt^ 327r4(l - e-'"?) , , 



47r^ + a^q^ 



aq a^g^(47r2 -f- a?q'^) 



where a is the width of the well. We thus obtain the 
dielectric functions e{q, d) as defined in Eqn. [31 where 
01 (q) remains as defined in Eqn. [2l giving: 



and where the ground plane is a 2D layer, we have: 



1 



e2D,ns,xf (97 d) 



= 1 



1-T2[1-G2((?)] 



(18) 



P2"'(q) - X2(q)[C/(q)pr*(q) + t/(q)pr'^(q)] (16) 

It is also important to account for the finite thick- 
ness of the screening and transport layers, which are 



J 



1 ^ (l-e-^'?'^)[l + TiGi(g)] 

emctai,s,xf (g, d) 1 - Ti(l - Gi{q) - e-29rf) 



(19) 



e2D,s,xf(g,ti) 



1 + TiGi(g) - T2[l - G2(g) - e-'^l + TiT2[Gi(g)G2(g) - Gi(g)(l ~ e'^^'')] 
[1 - Ti(l - Gi{q))][l - T2(l - G2{q))] - TiT2e-29'i 



(20) 



r 



where Ti = ')<^{q)F{q)V {q) and the additional subscript 
xf indicates the inclusion of exchange and finite thick- 
ness effects. As a consistency check, if we take the 2D 
screening layer to the metallic limit, by using G2{q) = 0, 
T2 = —q^^ I q and the limit q^^ — > 00, and return to zero 
thickness i^(g) = 1, then e2D,ns,xf in Eqn. [TBI reduces to 
em,etai,ns in Eqn. El and e2D,s,xf in Eqn. [201 reduces to 
tmetai,s in Eqn. 1111 If we separate the two layers by tak- 
ing d ^ 00 then both t2D,s,xj in Eqn.[20land emetaLs,xf 
in Eqn. ll9l rcduces to ^single in Eqn. 1131 We compare this 
work with related studies by Zheng and MacDonald^ in 
Appendix B. 



III. RESULTS AND DISCUSSION 

In this Section, we will use the various dielectric func- 
tions derived in Sect. II to answer an important physical 
question regarding recent experiments on screening long- 
range Coulomb interactions in 2D systems: Docs a 2D 
system screen as effectively as a metal when used as a 



ground-plane? 

We will answer this question in three stages. First 
we will consider the simplest possible case where there 
is no intralayer screening in the transport layer and the 
Thomas- Fermi approximation holds. Our results at this 
stage are directly applicable to ground-layer screening 
studies of dilute 2D systems, such as those investigating 
Wigner crystallization on liquid heliu m^^i^^ and the 2D 
insulating state in an AlGaAs/GaAs heterostructurCfi^ 
They may also be relevant to recent studies of the metal- 
insulator transition in Si MOSFETs^ where the gate 
is likely to produce significant ground-plane screening 
in the nearby 2DES located < 40 nm away, for exam- 
ple. Second, wc will then look at what happens when 
intralayer screening is introduced to the transport layer. 
This will allow us to understand why the ground-plane 
has such a significant effect on the insulating state in the 
experiment by Huang et ali^ and such little effect on the 
metallic state in the experiment by Ho et alJ^ Finally, 
since the experiments in RefsJ^ andi^ were performed 
at >> 1, we will investigate how our results change if 
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TABLE I: The d values for holes and electrons corresponding 
to the four cf"^ values considered in Sections III A and B. 





56.1 


9.89 


3 


1 


dholes (nm) 


50 


8.80 


2.67 


0.89 


^electrons (riHl) 


283.87 


50 


15.18 


5.06 



we extend beyond the Thomas-Fermi approximation and 
begin to account for finite thickness and exchange effects. 



A. Thomas-Fermi approximation in the absence of 
intralayer screening 

To get an understanding of the basic physics of our 
ground-plane screening model, we will begin by ignor- 
ing any effects of intralayer screening in the transport 
layer and use the Thomas-Fermi approximation to ob- 
tain the polarizability x(ci)- There are two important 
parameters in our equations: the layer separation d and 
the wave-number q, and to simplify our analysis we will 
make these parameters dimensionless by using q/q^^ and 
(P"^ = d X q"^^ hereafter. The Thomas-Fermi wave- 
number q'^^ contains all of the relevant materials param- 
eters involved in the experiment. In Ref. 14, where mea- 
surements were performed using holes in GaAs, = 12.8 
and m* = 0.38me, giving = 1.12 x lO^m^^ (i.e., 
io^^)^^ = 0.89 nm). The corresponding values for elec- 
trons with m* = 0.067me are q"^^ = 1.97 x lO^m"^ and 
(q^^)^^ = 5.06 nm). Table 1 presents the d values cor- 
responding to the four d^^ values that we will discuss in 
Sects. IIIA/B. The first two values correspond to d = 50 
nm for holes and electrons, the remaining two allow us 
to demonstrate what happens as the screening layer gets 
much closer to the transport layer in both cases. 

To facilitate a comparative analysis of the effective- 
ness of the screening, in Fig. [H^a) we plot the inverse di- 
electric function vs q/q^^ , with both a metal (solid 
blue/dotted green lines) and a 2D system (dashed red 
lines) as the screening layer, for the four different d"^^ 
values listed in Table 1. Note that for the metal screening 
layer case {emetai.,nsiq,d) in Eqn. [5]), we have explicitly 
parameterized q and d into q/q^^ and , in order to 
plot the metal and 2D screening layer cases on the same 
axes. The metal data for d^^ = 56.1 is presented as 
dotted green line as it serves as reference data for later 
figures. Note that e"^ ~ \ corresponds to no screening, 
and = corresponds to complete screening of a test 
charge placed in the transport layer. Considering the 
large d^^ limit first, e^^ only deviates from 1 at small 
qjq^^ 1 and heads towards e^^ = as qjq^^ — > 0. In 
other words, screening is only effective at large distances 
from a test charge added to the transport layer. This 
makes physical sense if one considers the electrostatics of 
ground-plane screening. The ground-plane acts by inter- 
cepting the field lines of the test charge such that they are 
no longer felt in other parts of the transport layer. This 
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FIG. 3; 

(a) The inverse dielectric function e~^(g,d) vs q/qTF with 

no intralayer screening in the transport layer. Data is 
presented for metal (solid blue and dotted green lines) and 
2D (dashed red lines) screening layers for the four dxF 
values presented in Table 1. The metal gate data for 
dxF = 56.1 appears as a dotted green line as it serves as 
reference data for later figures, (b) The relative effectiveness 
of the ground-plane screening due to a 2D screening layer 
compared to a metal screening layer, as quantified by the 
ratio Cetai.ns/f^D.ns. ^s q/qTF for the four drF values. 



is only effective at distances from the test charge that 
are much greater than the ground-plane separation d, 
and thus the ground-plane acts to limit the range of the 
Coulomb interaction in the transport layer, as pointed 
out by Peeters.— With this in mind, it is thus clear why 
the point of deviation from = 1 shifts to higher values 
of q/q^^ as d^^ is reduced. Indeed, all four lines pass 
through a common value when q/q^^ = jt-t reflect- 



ing this electrostatic aspect of ground-plane screening. 

Turning to the central question of the effectiveness of 
a 2D layer as a ground-plane, in Fig. EJa) it is clear 
from the increasing discrepancy between the solid and 
dashed lines that the 2D system becomes less effective 
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than a metal as d?"^ is reduced. To quantify this, in 
Fig. El^b) we plot the ratio of the two dielectric constants 
^metai ns/^2D ns' with a ratio of 1 indicating equivalent 
screening and < 1 indicating that a 2D system is less ef- 
fective than a metal. For large separations, for example 
d^^ = 56.1, which corresponds directly to the experi- 
ment by Ho et ai, a 2D system screens as effectively as 
the metal gate to within 1%. However if the screening 
layer is brought very close to the transport layer cf^^ ~ 1 
(i.e., the screening layer is only a Thomas-Fermi screen- 
ing length away from the transport layer) then the effec- 
tiveness of the 2D system as a ground-plane is reduced to 
~ 66% of that of a metal layer at an equivalent distance. 
It is important to note that correlations between the two 
layers can be significant for such small separations, and 
hence this increasing discrepancy should be considered as 
a qualitative result only. Furthermore, as we will sec in 
Sect. niC, exchange actually acts to enhance the effec- 
tiveness of the 2D system as a ground-plane, making the 
Thomas-Fermi result above a significant underestimate 
of the true ground-plane screening of a 2D system in the 
low density limit. 

B. Thomas- Fermi Approximation with Intralayer 
Screening in the Transport Layer 

We now add intralayer screening in the transport layer 
to our Thomas-Fermi model, and begin by asking: What 
is the magnitude of this intralayer screening contribu- 
tion, independent of any ground-plane screening effects? 
In Fig.mja), we plot the inverse dielectric function CgjjjgiQ 
(dash-dotted black line) for a 2D system with intralayer 
screening and no nearby ground-plane. For comparison, 
we also show the data from Fig. for a metal ground- 
plane with d-TF = 56.1 (green dotted line) and the ex- 
pectation with no screening, e^^ = 1 for all q/qxp (grey 
dashed horizontal line) in Fig. HJa). It is clear that the 
addition of intralayer screening has a very significant im- 
pact on the dielectric function, more so than the addition 
of a ground-plane. Indeed, returning to an electrostatic 
picture and ignoring exchange and correlation effects, 
^single should assume the (f^ — > limit of e^D^g, the 
Thomas-Fcrmi model in the absence of intralayer screen- 
ing. 

We now reintroduce the ground-plane, and in 
Fig. m (a), we plot the combined screening contribu- 
tions for metal i^^t^^i si<l,d), solid blue lines) and 2D 
(e^Q j,((7, d), solid red hues) screening layers. These are 
shown for the four different values of cP"^ in Table 1. 
The values for the 2D system are offset vertically by 
—0.2 for clarity. The intralayer screening and ground- 
plane screening both contribute to the total screening, 
albeit on different length scales. This can be seen by 
comparing the data in Fig. [UJa) to that in Fig. [3l^a), 
with the intralayer screening clearly the dominant con- 
tribution. As a result, distinguishing between individ- 
ual traces in the sets corresponding to the metal (blue 




^TF^ 5g .| ^mftai(^''^)' metal ground plane 

— — RopCq.d), 2D ground plane 
J I I I I I I I I I 

0.001 0.01 0.1 1 



FIG. 4: 

Effect of a screening layer on a 2D system with intralayer 
screening, (a) Firstly, in order to show the relative effects of 
intralayer and ground plane screening, we plot the dielectric 
functions = 1 corresponding to no intralayer or ground 
plane screening (grey dashed horizontal line) , ^^j, with 

metal screening layer at (Itf = 56.1 and no intralayer 
screening (dotted green line, data from Fig. |3ja)), and e^Jj^gj^ 
with intralayer screening but no ground-plane (dash-dotted 

black line). We then consider the effect of the metal 
screening layer when the intralayer screening is included, by 
plotting g(q,d) (solid blue lines) for the four values of 

d^^ shown in Table I. Moving through the traces from 
upper left to lower right corresponds to decreasing d^^ . A 
similar set of curves are shown for the case of a 2D screening 

layer (e^^ ^{q,d), solid red lines) which has been offset 
vertically by -0.2 for clarity, along with a duplicate of e^^gi^ 

(dash-dotted black line). Since the dielectric functions 
almost lie on top of each other when intralayer screening is 
present, in (b) we plot R, the relative enhancement of 
due to the ground plane. Calculations for the four different 
values of drp in Table 1 are shown for metal (solid blue 
lines) and 2D (dashed red lines) ground-planes. 
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lines) or 2D (red lines) ground planes in Fig. lUJa) is dif- 
ficult. Hence, to better quantify the enhancement that 
the ground-plane gives over intralayer screening alone, in 
Fig. [IJb) we plot the relative ground-plane enhancement 

ii'mctaU = (e.4nglc-emetal,s)/|es7nglel (''OM bluC huCs) and 

^2D,s = (esi,igie - em,s)/|es"Ingiel (dashed red lines), re- 
spectively. Note that the ground-plane only provides sig- 
nificant enhancement over intralayer screening alone as 
becomes small, and as in Sect. IIA, only provides 
enhancements at small q/q'^^. The small discrepancies 
between the data for the metal and 2D screening layers 
in Figs. IDJb) directly reflect the increased effectiveness of 
the metal ground-plane over a 2D ground-plane shown in 
Fig.Hb). 

The data in Figs.[3]and [Uprovide an interesting insight 
into recent experiments on ground-plane screening in 2D 
hole systems in the insulating and metallic regimesJ^ii^ 
Due to the low hole density and conductivity in the insu- 
lating regime, intralayer screening is less effective and the 
dominant contribution to screening is the ground-plane, 
which acts to limit the length scale of the Coulomb inter- 
actions, as Fig. Efa) shows. This results in the ground- 
plane having a marked effect on the transport properties 
of the 2D system, as shown by Huang et alJ^ In com- 
parison, for the metallic state, where the density and 
conductivity are much higher, intralayer screening is the 
dominant contribution, and a ground-plane only acts as 
a long-range perturbation to the screening, as shown in 
Fig. BJa). This perturbation to the intralayer screening 
is particularly small at (f"^ = 56.1 and results in the 
ground-plane having relatively little effect on the trans- 
port properties in the metallic regime, as found by Ho et 
alJ^ Although Fig. |4] (b) suggests that decreasing cf^ 
will increase the effect of the ground plane, in practice 
there are issues in achieving this. For holes in GaAs, 
there is little scope for further reducing d due to increas- 
ing Coulomb drag and interlayer tunnelling effects. Also, 
in our model we have neglected interlayer exchange and 
correlation effects, and these may become significant at 
these lower distances. 



C. Beyond the Thomas-Fermi Approximation 
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FIG. 5; 

(a) Plots of vs q/qxp for a metal (solid blue and dotted 
green line) and 2D (dashed red lines) ground plane for 
d = 30 and 50 nm accounting for exchange and finite 
thickness effects but ignoring intralayer screening in the 
transport layer. The dotted green line corresponds to that in 
Fig. 131 a), (b) A plot of the relative screening effect of a 2D 
layer compared to a metal (solid lines), as quantified by the 
ratio eZ}t Jf-'^-^- In contrast to the results for the 

metal/ Ilj 

Thomas- Fermi model (dashed line - data from Fig. [3jb)), we 
find that a 2D layer is actually more effective than a metal 
as a ground-plane when exchange and finite thickness effects 
are included in the calculation. 



Following our relatively simple treatment of ground- 
plane screening above, it is now interesting to ask how the 
results of our calculations change if we extend our model 
to account for two phenomena ignored in our Thomas- 
Fermi model: exchange effects at low densities, and the 
finite thickness of the screening and transport layers. 

The inclusion of the Hubbard local field correction 
G{q), finite thickness form-factor F{q), and the use of the 
Lindhard function for x^W) adds two new parameters to 
the analysis, the well thickness a = 20 nm and the Fermi 
wave- vector kp. This removes our ability to reduce the 
problem down to a single adjustable parameter dP"^ as 
we did in Sects. IIIA and B. Additionally, accounting for 
finite well width puts a lower limit on rf, which must be 



greater than a to ensure that the wells remain separate. 
Hence for the remaining analysis we will only consider 
d = 50 and 30 nm, which correspond to (f"^ = 56.1 
and 33.7, respectively. As in earlier sections, we will first 
analyse the dielectric function ignoring intralayer screen- 
ing in the transport layer, which we achieve by setting 
pf'^ = 0. 

In Fig. [Slja) , we plot ejj] ,jg (dashed red lines) ob- 
tained using Eqn. [201 for d = 50 and 30 nm, and for 
comparison, e~otai ns ^'^^ d = 50 nm (dotted green line) 
from Fig. El^a) and the corresponding result for d = 30 
nm (solid blue line). One of the more significant effects 
of exchange in 2D systems is that it leads to negative 
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compressibilit}i2^i2i for > 2. It is well known from the 
field penetration experiments of Eisenstein et al^ that 
negative compressibility is related to an overscreening of 
the applied electric field, leading to a negative penetra- 
tion field. In our calculations, a similar overscreening 
is observed, with a 2D system producing more effective 
ground-plane screening than a metal gate at intermedi- 
ate qjq^^ in Fig. O^a). The enhanced screening when 
the ground-plane is a 2D system is evident in Fig. [5jb), 
where we plot the ratio e"^^^^^ ns/^2D ns' which takes val- 
ues greater than 1 for q/q^^ ^ O-l- 

We now reintroduce intralaycr screening in the trans- 
port layer, and in Fig. [B]^a) we plot e.'^^^^^gx g (solid blue 

line) and ejj] ^ (dashed red line) for d = 50 and 30 nm. 
The values for the 2D system arc offset vertically by —0.2 
for clarity. For comparison, we also plot e^ngi,, (dash- 
dotted black lines — duplicated and offset vertically by 
—0.2), along with the data from Fig. [Slja) for a metal 
ground-plane at d = 50 nm with no intralayer screening 
(green dashed line), and the expectation with no screen- 
ing ~ 1 for all q/qrp (grey dashed horizontal line). 
As wc found earlier with the Thomas-Fermi model (see 
Fig.Ua)), the inclusion of intralayer screening has a pro- 
found effect on the dielectric function, contributing sig- 
nificantly more to the overall screening than the addition 
of a ground-plane does alone. This demonstrates the ro- 
bustness of one of the key results of Sect. IIIB, namely 
that in the metallic regime,— where intralayer screening 
effects are significant, the ground-plane screening contri- 
bution is overwhelmed by the intralayer screening con- 
tribution. This leads to a significantly reduced ground- 
plane effect than one would expect from studies in the 
insulating regimeii^ 

The effect of including exchange and finite thickness 
effects in the calculation is evident by comparing ^ 

in Fig. IHa) with g in Fig. SJa). Considering the 
individual contributions, because F{q) < 1, the finite 
thickness of the quantum well acts to reduce the effec- 
tiveness of the 2D layer as a ground-plane. In contrast, 
the negative compressibility produced by the exchange 
contribution acts to significantly enhance the screening, 
and as Fig. [DJa) shows, has its most significant impact 
at intermediate q/qxp, where the dielectric function be- 
comes negative, as discussed by Dolgov, Kirzhnits and 
MaksimoA^, Ichimaru^i, and Iwamoto^. The combined 
effect of G{q) and F{q) is to significantly enhance the 
screening at intermediate q/qrp whilst reducing it to 
levels comparable to the metal ground-plane for large 
q/qTF- III other words, the added density-dependence in 
our Hubbard model leads to enhanced mid-range screen- 
ing at the expense of short-range screening. A physi- 
cal interpretation for this behavior is that at low den- 
sities there are insufficient carriers available to screen 
effectively close to a test charge, whilst at intermedi- 
ate ranges, the negative compressibility produced by ex- 
change leads to a higher availability of carriers and better 
screening than there would otherwise be at higher car- 
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FIG. 6: 

Effect of a screening layer on a 2D system with intralayer 

screening, with exchange and finite thickness efi'ects 
included, (a) Firstly, in order to show the relative efi'ects of 
intralayer and ground plane screening, we plot the dielectric 
functions = 1 corresponding to no intralayer or ground 
plane screening (grey dashed horizontal line) , e~etai ns with 
metal screening layer at d = 50nm and no intralayer 
screening (dotted green line, data from Fig. |3ja)), and 
^single xf with intralayer screening but no ground-plane 
(dash-dotted black line). We then consider the efltect of the 
metal screening layer when the intralayer screening is 
included, by plotting ^ ^{{q, d) (solid blue lines) for 

d = 30 and 50nm. A similar set of curves are shown for the 
case of a 2D screening layer (e^^ ^ xf '^)' solid red lines) 
which has been ofiset vertically by -0.2 for clarity, along 
with a duplicate of e^ngj^ (dash-dotted black line) . Since 
the dielectric functions almost lie on top of each other when 
intralayer screening is present, in (b) we plot R, the relative 
enhancement of due to the ground plane. Calculations 
for d = 30 and 50nm are shown for metal (solid blue lines) 
and 2D (dashed red lines) ground-planes. 

rier densities where exchange is not as significant. It is 
also interesting to consider why the introduction of ex- 
change and finite thickness effects have such a profound 
effect on the intralayer screening contribution compared 
to the ground-plane screening contribution. This occurs 
because the impact of G{q) and F{q) on the ground-plane 
contribution is strongly attenuated by the e"^*"^ terms 
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that appear in Egns. [T8l and [20l Such terms don't occur 
for the intralayer screening contribution, which signifi- 
cantly enhances the impact of the negative compressibil- 
ity, as is clear by comparing Fig. [5Ja) with Fig. Efa). 

We close by considering the relative effectiveness of 
the metal and 2D ground-planes with all considera- 
tions included in the calculations. In Fig. ^h) we plot 
the relative ground-plane enhancements i?mctai,s,xt = 
(e.^^gicxf - emctai,s,xf)/|es"Iiigio,xfl (soHd blue lines) and 
i?2D,s,xf - (Csingicxf - e^D,s,xf)/k.^^gic,xfl (dashcd red 
lines) for d = 50 and 30 nm. As in Fig.[5t^b), we find that 
exchange, finite thickness and intralayer screening result 
in the 2D ground-plane screening more effectively than 
a metal ground-plane, with the difference between the 
two becoming greater as d is decreased. For = 50 nm, 
the ground-plane separation used in Ref. 14. the ground- 
plane has significantly more effect (~ 8 — 9%) than it 
does in the more simple Thomas-Fermi model 1%) 
presented earlier. 

IV. SUMMARY OF RESULTS AND 
COMPARISON WITH EXPERIMENT 

We have performed theoretical calculations to investi- 
gate the relative effectiveness of using a metal layer and 
a 2D system as a ground-plane to screen Coulomb in- 
teractions in an adjacent 2D system. This is done for 
two cases: the first is the relatively simple Thomas-Fermi 
approximation, and the second is the Hubbard approx- 
imation, where we account for exchange and also finite 
thickness effects. This study was motivated by recent 
experiments of the effect of ground-plane screening on 
transport in semiconductor-based 2D systems. 

There were three key findings to our study. Firstly, 
a 2D system is effective as a ground-plane for screening 
Coulomb interactions in a nearby 2D system, which was 
an open question following the recent experiment by Ho 
et ali^ In the Thomas-Fermi approximation, a metal and 
a 2D system are almost equally effective at screening the 
long-range Coulomb interactions in the nearby 2D sys- 
tem, with the metal becoming relatively more effective 
as the ground-plane separation d is decreased. 

Secondly, our calculations provide an explanation for 
why ground-plane screening has much more effect in the 
insulating regime than it did in the metallic regime. Due 
to the low hole density and conductivity in the insulat- 
ing regime, intralayer screening is weak and the domi- 
nant contribution to screening is the ground-plane, which 
acts to limit the length scale of the Coulomb interactions. 
This results in the ground-plane having a marked effect 
on the transport properties of the 2D system, as shown by 
Huang et alr^ In the metallic regime, intralayer screen- 
ing cannot be ignored. In addition to being the dominant 
contribution for long-range interactions (i.e., at small q), 
the intralayer screening contribution is non-zero over a 
much wider range of q, turning the ground-plane contri- 
bution into little more than a small change to the overall 



screening in the 2D system, which is consistent with the 
experiment by Ho et alr^. 

Finally, since both experiments were performed at 
Ts » 1, where the Thomas- Fermi approximation is in- 
valid, we reconsider our calculations involving 2D sys- 
tems using the Hubbard approximation for the local field 
correction. We show that our argument regarding the 
physics of ground-plane screening in the metallic and 
insulating regimes remains robust, but that exchange 
effects lead to a 2D system being more effective than 
a metal layer as a ground-plane. This is due to the 
exchange-driven negative compressibility that occurs22, at 
rs>2. 



V. FURTHER WORK 

While our results suggest that ground-plane effects 
on a metallic transport layer should strengthen as the 
ground-plane separation d^^ is reduced, there are a num- 
ber of issues that complicate this argument. Firstly, for 
holes in GaAs, such as the experiment in Ref. 14, there 
is little scope to further reduce d due to the increas- 
ing Coulomb drag and interlayer tunnelling effects that 
would result. However, it may be possible to experi- 
mentally modify d^^ by moving to a different material 
system where {q'^^)~^ is larger. For example, in InAs^, 
where m* = 0.026TOe and e = 14.6eo, we would have 
(q^P)-^ 14.9 nm, or InSb^ where m* = 0.0145™^ 
and e = 17.7eo gives (g^^)^^ = 32.3 nm. These (q^^)'^ 
values are 17 and 28 times larger than those in Ref. 14, 
respectively. This would allow us to reduce dxp without 
changing d, thus avoiding the problems above. 

We note that our model neglects interlayer exchange 
and correlation effects, which may become significant at 
these small distances d, as suggested by calculations at 
Ts = 4 by Liu et al. — It could be interesting to inves- 
tigate the effect of including the interlayer exchange and 
correlation effects on the ground plane screening at small 
d. This may require using better approximations for the 
local field correction such as that developed by Singwi, 
Tosi, Land and Sjolander^^ (STLS), as there is no equiv- 
alent to the Hubbard approximation for interlayer local 
field corrections. 

We also note that using the technique in Ref. 14 and 
the theory presented here, it would be possible to study 
the breakdown of intralayer screening in the transport 
layer as it is evolved from the metallic to insulating 
regime. This could be compared with compressibility 
measurements of a 2D system across the apparent metal- 
insulator transition^, possibly providing new insight into 
the mechanism driving this transition. 

Lastly, in this paper we only calculate the screening of 
the ground-plane on the transport layer via the dielectric 
function. It would be interesting to take this work fur- 
ther to calculate the effect of the ground-plane on the ac- 
tual carrier transport through the transport layer. Com- 
bining the theory presented here, and various models of 
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the metallic and insulating behaviours (see review pa- 
pers^iii^), it may be possible to determine how each of 
the models are affected by the presence of a ground plane, 
and would allow us compare this with the experimental 
data in more detail. 
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APPENDIX A: BRIEF REVIEW OF SCREENING 
THEORY FOR A SINGLE 2D SYSTEM 

In this Section we briefly review the basics of screening 
in a single 2D system. Readers familiar with screening 
theory may wish to proceed directly to Sect. II. A more 
extended discussion can be found in Refs.— i^^ii^. 

Screening occurs when the carriers in a 2D system re- 
organize themselves in response to some added 'exter- 
nal' positive charge density, leading to an electrostatic 
potential determined by Poisson's equation. This reor- 
ganization produces a negative 'induced' charge density 
that acts to reduce or 'screen' the electric field of the 
external charge. In proceeding, it is mathematically con- 
venient to instead treat the problem in terms of wave- 
vectors (g-space) so that the (intralayer) Coulomb po- 
tential V{r) = becomes V{q) = 

There are two key parameters of interest in an analysis 
of screening. The first is the polarizability x{q)7 which 
relates the induced (screening) charge density p™''(q) to 
the external (unscreened) potential (/i'^'^*(q): 



P"'(q)-X(g)0^"*(q) (Al) 

The second is the dielectric function e(g), which relates 
the total (screened) potential 0(q) to the external (un- 
screened) potential (/)'^'^*(q): 



0^-*(q)/e(g) 



(A2) 



Conceptually, the polarizability describes how much in- 
duced charge density is produced in response to the addi- 
tion of the external charge density, hence it is also often 
called the density-density response function;^ The di- 
electric function is a measure of how effective the screen- 
ing is: = 1 corresponds to no screening and = 
corresponds to perfect screening.— The two parameters 
can be linked via (j)'^^* and the Coulomb potential V{q), 
such that: 



1 



— = 1 + V{q)x{q) 



(A3) 



The results above are precise aside from the assumption 
of linear response. However, continuing further requires 



calculation of x{q)- This cannot be achieved exactly, and 
requires the use of approximations. In the simplest in- 
stances, a combination of the Thomas-Ferm i^^'^° (TF) 
and Random Phase Approximationsii (RPA) can be 
used. However, to properly account for exchange and/or 
correlation, particularly at lower carrier densities, more 
sophisticated approximations, such as those developed by 
Hubbard2i or Singwi, Tosi, Land and Sjolanderi^ (STLS) 
should be used. For a single 2D layer, this leads to a cor- 
rection to the induced charge: 



P"'(q) = X°iqW\^) + V{q)r''{q){l - G{q))] (A4) 
where G{q) is the local field factor. This results in: 



X"iq) 



1-V{q)x'^{q)[l-G{q)] 



(A5) 



The local field factor can be calculated in nu- 
merous ways.— In this work, we use the Hubbard 
approximation2Zi2^, which gives a local field factor: 



G{q) 



2vV 



(A6) 



where kp = \/27rp is the Fermi wave- vector. Al- 
though better approximations are available^, the Hub- 
bard approximation is sufficient to introduce a density- 
dependence into the screening, unlike the Thomas-Fcrmi 
approximation, which is density- independent. 



APPENDIX B: COMPARISON WITH OTHER 
WORK ON BILAYER SCREENING 

In this Appendix, wc discuss how the analytical ex- 
pression we obtain for ^ xf compares with other 
works on linear screening theory for bilayer 2D systems 
produced in double quantum well heterostructures, in 
particular, that of Zheng and MacDonald^i. Note that 
we have translated the equations from Ref. 21 into the 
notation used in our paper for this Appendix. 

Zheng and MacDonald begin by defining a density- 
density response function (polarizability) Xij "w^) for 
their bilayer 2D system by: 



p,{q,Lu) = ^x^,('7,'^)0r('?''^) (Bl) 
j 

where p is the linear density response (i.e., induced charge 
density), 4>j^* is the external potential, and i,j = 1,2 
are the layer indices with 1 being the transport layer 
and 2 being the screening layer. Zheng and MacDon- 
ald then use the random phase approximation^^ (RPA) 
and Singwi, Tosi, Land and Sjolandcr (STLS) approxi- 
mation^ to obtain an expression for the polarizability: 
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U{q)[G2iiq)~-l] 



r 



f/(9)[Gi2(g)-i] 

[x°(<7,^)]-i-n9)[l-G22(<?)] 



(B2) 



where Gij{q) are the local field factors that account for 
the effects of exchange and correlationi^ For compari- 
son with our work, we will consider w = 0, and ignore 
interlayer exchange and correlations by setting (g) = 
G2i(g) = 0, Gii(g) = Gi{q) and G22(g) = G2(g). The 
latter approximation will be valid for large d, but we 
would expect that Gij would become more significant at 
lower distances. This is seen in the work of Liu et al. 
—, in which Gu and Gy are calculated using STLS for 
different d, at = 4. 

In our work, we are seeking to obtain an effective sin- 
gle layer dielectric function for the transport layer only. 



J 



Hence we only put external charge density pf^* (q) in the 
transport layer and set the external charge density in the 
screening layer P2^*((?) to zero. This results in external 



potentials in the two layers of (/if (q) = V{q)pi{q), and 



lext { 
"2 



'(g) = U{q)pf^*{q). The total potential in the trans- 
port layer can thus be expressed as: 



Ml) = + V{q)pT{q) + U{q)pT{q) (B3) 

For the dielectric function of the transport layer, as de- 
fined in Eqn. [31 this results in: 



e{q,d) 



1 + V{q)xil{q) + U{q)xi2{q) + U{q)x2l{q) + e^'^'^U {q)x22{q) 



(B4) 



This is analogous to Eqn. lA3l for the single layer case. In- 
deed, by applying li ^ oo to Eqn. lB4l reduces to Egn. lASl 
Finally, obtaining the matrix elements Xij [q) by inverting 
Eqn. IB2I and inserting them into Eqn. IB41 we obtain the 
same expression as that given for e^u sxfil) ™ Eqn. [^D] 
after returning to zero thickness (i.e., F{q) = 1). 



APPENDIX C: COMPARISON WITH THE 
SCREENING OF A PERPENDICULAR 
ELECTRIC FIELD BY A 2D SYSTEM 

In this Section we show that the calculations in this 
paper, which describe the screening of in-plane point 
charges (and also arbitrary in-plane charge distributions) 
by a 2D system and also an adjacent ground plane, are 
consistent with the equations by Eisenstein et al^ de- 
scribing the penetration of a perpendicular electric field 
across a 2D system. 

In order to calculate the penetration of the perpendic- 
ular electric field, it is necessary to consider a slightly 
different configuration than previously used in Fig. (Hb). 
Figure [7] shows a schematic of the system we now con- 
sider. We still have two 2D systems (labelled 1 and 2) 
separated by a distance d. We now have no external 
charge in either of these layers, and only induced charges 
p™'' and p™''. In order to apply an electric field across 
layer 1, we place a layer of external charge p^* a distance 
d2 below layer 1. The net electric fields Eq and Ep are 
shown in Fig. 7, and have been related by Eisenstein et 
al. using Eqns. (5) and (6) in Ref. 34. Translating into 



2D Screening layer 2 { 02, p") 



t t t E, 



2D Screening layer 1 { <X>,, p,'"^] 
A k k k 



f + + + 4 



+ 



2D layer {^o, /?o"'} 



FIG. 7; 

Schematic showing the configuration used to calculate the 
penetration field across a 2D system. 



the notation of our paper, and using = ^e^(f^)i, this 
results in: 



dE,, 



dEo -e(x'i' + x'i) + dxU'i 



(CI) 



We now try to calculate the electric fields Eq and Ep 
directly using our model. The induced charges p™** a-nd 
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Pj"'^ are calculated in a similar fashion to previously in 
Section II. C, although Eqns. [15] and [H] need to be 
modified slightly to take into account the position of the 
external charge. This results in: 



pr'(q) = xi(q)[0r*(q) + t/(q)pr'(q)] 
pr'(q) = X2(q)[0r*(q) + t^(q)pr'(q)] 



where 



lext 



(q) 



F(q)e-i'^Vo"*(q) 



q) = V{q)e- 



-C{{d-\-d2} ^GXt 



Po"*(q)- 



(C2) 
(C3) 



(C4) 
(C5) 



We set pg^*(r) = co, where o-q is a constant 2D sur- 
face charge density. This simplifies our model into the 
ID problem considered by Eisenstein et al^ In g-space 
this gives pg'^*(q) = CTo(27r)^(5(q), where S{q) is the Dirac 
delta function. Solving equations [U2] and [CS] simultane- 
ously we obtain p*"''(q) — cri(27r)^(5(q) with: 



CTi = 



(72 



CToeX2 



xl) 



dxlx'l 



(C6) 
(C7) 



where we have used the property of delta functions that 
/(q)(5(q) = /(0)(5(q) for an arbitrary function /(q). We 
note that the above equations IC6I and IC7I are valid even 



for an arbitrary local field correction, as at g = the 
Lindhardt function is equal to the Thomas-Fermi func- 
tion, and for any self-consistent local field correction 
G(0) = (see equation 3.25 in Ref. 37.) Similarly, the 
eauations lC6l and lC7l are valid even when allowing for the 
finite thickness of the 2D systems, since for an arbitrary 
form factor F{q) we have F(0) = 1. 

We can now calculate the electric fields Eq and Ep. 
For a single layer of uniform 2D charge density cr, the 
perpendicular electric field is given hy E = By con- 
sidering the positions of the electric field Eq and Ep with 
respect to the charge layers pq, pi, and p2, we obtain: 



En 



En = 



which results in 



E^, 



p 

Eq 



-I- (Tl — 0"2 

2e 

gp - O-j - (7-2 

2e 



-ex'2 



(C8) 
(C9) 



(CIO) 



We note that in this model, the ratio is equal 

to as we have neglected other charges commonly 

present in 2D systems, such as regions of modulation 
doping. Our resulting equation lClOl is thus in agreement 
with the equations used by Eisenstein et al^ (equation 

EH. 
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